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The discovery and timing of radio pulsars within the Galactic centre is a fundamental aspect of 
the SKA Science Case, responding to the topic of “ Strong Field Tests of Gravity with Pulsars 
and Black Holes ” (Kramer et al. 2004; Cordes et al. 2004). Pulsars have in many ways proven 
to be excellent tools for testing the General theory of Relativity and alternative gravity theories 
(see Wex (2014) for a recent review). Timing a pulsar in orbit around a companion, provides a 
unique way of probing the relativistic dynamics and spacetime of such a system. The strictest 
tests of gravity, in strong field conditions, are expected to come from a pulsar orbiting a black 
hole. In this sense, a pulsar in a close orbit (P OY b < 1 yr) around our nearest supermassive black 
hole candidate, Sagittarius A* - at a distance of ~ 8.3 kpc in the Galactic centre (Gillessen et al. 
2009a) - would be the ideal tool. Given the size of the orbit and the relativistic effects associated 
with it, even a slowly spinning pulsar would allow the black hole spacetime to be explored in great 
detail (Liu et al. 2012). For example, measurement of the frame dragging caused by the rotation 
of the supermassive black hole, would allow a test of the “cosmic censorship conjecture.” The 
“no-hair theorem” can be tested by measuring the quadrupole moment of the black hole. These 
are two of the prime examples for the fundamental studies of gravity one could do with a pulsar 
around Sagittarius A*. As will be shown here, SKA1-MID and ultimately the SKA will provide 
the opportunity to begin to find and time the pulsars in this extreme environment. 
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1. Background 

A combination of radio and infrared (IR) astrometric observations has established the pres¬ 
ence of a dark mass in the Galactic centre (GC), whose properties are most consistent with being a 
supermassive black hole (BH) (Eckart & Genzel 1996; Ghez et al. 1998; Reid & Brunthaler 2004; 
Ghez et al. 2008; Gillessen et al. 2009a,b). Known by its radio nomenclature of Sagittarius A* 
(Sgr A*), current estimates are that the mass of this BH is 4.30 ± 0.20( stat ) ±0.30( sys ) x 10 6 M 0 
(Gillessen et al. 2009a). Monitoring the orbits of the closest stars SO-2 (Schodel et al. 2002) and 
S0-102 (Meyer et al. 2012) could provide constraints on the extent to which their - orbits are per¬ 
turbed which in turn would constrain the enclosed mass, which could be stars, dark matter, or both 
(Sabha et al. 2012; Gillessen et al. 2009b). Current work with the Keck Telescope, the VLT and 
with the VLT interferometer (VLTI) and proposed work with extremely large telescopes (ELTs; 
e.g. the European ELT and the GSMT) will lead to IR astrometric precisions of 0.5 mas with single 
apertures and 10 /./as with the VLTI and the GRAVITY instrument (Weinberg et al. 2005; Eisen- 
hauer et al. 2011). These resolutions allow orbital perturbations of 4.15 and 0.08 au, respectively, 
to be measured that will arise from General Relativity (GR) and from Newtonian effects caused 
by other stars and dark matter. It is expected that source confusion in IR surveys with ELTs will 
limit the sample of stars to those with orbits larger than about 200 to 300 au; ~ 100 stars with 
semi-major axes < 3000 au for a 30-m ELT (Weinberg et al. 2005). 

The discovery and timing of radio pulsars orbiting Sgr A*, and within the central stellar clus¬ 
ter, will provide a powerful means to test GR, as well as probe the dark matter content and magne¬ 
toionic medium within the nucleus of our Galaxy. An SKA radio pulsar program can yield much 
more precise localizations on objects arbitrarily close to the innermost stable circular orbit (ISCO, 
radius >0.1 au, depending on the spin of the Sgr A* BH). Confusion limits do not apply to the 
time-varying signals of pulsars. In comparison with the best forecasted precisions of IR astrometry 
and Doppler measurements, pulsar timing has the potential of increasing the testing and diagnostic 
precision by at least three orders of magnitude. Pulsars offer this prospect because timing pre¬ 
cisions correspond to spatial localizations (from ranging ) smaller than a light second (0.2 /./as) 
compared to the localizations expected from the best proposed IR interferometry, of no better than 
10 jUas (0.1 au). 

2. The pulsar population in the GC 

Several lines of argument suggest that not only is there a large neutron star (NS) population in 
the GC but that many NSs should be active radio pulsars. Figure [I] (left panel) shows schematically 
the contents of the GC as determined both empirically and from theoretical arguments about galaxy 
formation. Ghez et al. (2008) and Gillessen et al. (2009b) show that the stellar cluster around 
Sgr A* consists mainly of early-type stars, whose short lifetimes (~ 10 7 yr) means that either 
massive star formation in the region is ongoing or that there has been a starburst of order 6 Myr 
ago, although the presence of young, massive stars in a region where tidal forces should inhibit star 
formation present a conundrum: “the paradox of youth” - (Ghez et al. 2003). Nonetheless, stars 
with masses of 10-20 M 0 are NS progenitors, and the presence of a population of young massive 
stars near Sgr A* suggests that there must exist a numerous NS population as well. The time scales 
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Figure 1: (Left) Schematic view of the matter content of the GC. Young stars (e.g., SO-2) are known to 
approach as close as 0.1 " (1000 AU) to Sgr A*. The combination of young, hot stars in the GC and mass seg¬ 
regation suggests that there are likely to be compact stellar remnants (NSs and BHs) at least within 1000 au, 
and potentially much closer. Finally, there could also be dark matter particles. The SKA, operating at fre¬ 
quencies near 10 GHz, will search for NS, in the form of active radio pulsars, as close as possible to Sgr A*. 
(Right) Depth of gravitational potential ( GM/ac 2 ) probed by an orbit of semi-major axis a as a function of 
the mass of the binary system M. The regions probed by solar system tests and known pulsars in binary 
systems are indicated, with the double neutron star systems PSR B1913+16 and PSR J0737—3039 shown 
explicitly. The vertical cyan line indicates the mass of Sgr A* (GC-BH). The red line segment shows the 
region probed by the S stars in the GC. It is clear that substantially deeper gravitational potentials (larger 
GM/ac 2 ) remains to be probed, and would be accessible to radio pulsar observations (Kramer et al. 2004). 


also work out for many of these NSs to be active radio pulsars: radio emission of canonical 10 12 G 
objects lasts for ~ 10 Myr, the time scale for spindown to lengthen spin periods past the “death 
line” of a few seconds. 

Estimates of the supernova rate within 100 pc of Sgr A* based on the number of young 
stars in the region and on the heating of X-ray gas range from 10 3 yr~* to 10 2 yr _1 and are 
consistent, after accounting for the fraction of supernovae that produce pulsars, with the number 
of point radio sources in the region that could be pulsars (Lazio & Cordes 2008). Such agreement 
is coarse but suggestive that there are of order 10 3 to 10 4 active pulsars in the region, of which 
about 200 to 2000 (20%) would be beamed toward the Earth. More recently, and by combining 
results from the GC radio point source counts described above, the detected GC pulsar population 
at that time, non-detections in high-frequency pulsar surveys of the inner parsec, radio and gamma- 
ray measurements of diffuse emission, infra-red observations of massive stellar GC populations, 
pulsar wind nebulae candidates, and the supernova rates from X-ray data, Wharton et al. (2012) 
estimate that ~ 10 3 active radio pulsars in the inner parsec should be beamed toward Earth. 

Pulsars in the GC, like other stellar populations, will orbit the central BH Sgr A* mostly as 
isolated objects but some will be in binaries with companions of all types: main-sequence and 
post-main-sequence stars, white dwarfs, NSs, or even BHs (Faucher-Giguere & Loeb 2011). A 
population of “recycled” pulsars (i.e., NSs spun up to millisecond periods by accretion) is also 
expected because the stellar density in the GC is large enough so that tidal capture and exchange 
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Figure 2: The positions of detected radio pulsars in the central 0.5° toward the GC overlaid on a 10.55 GHz 
continuum map made with the Effelsberg telescope (Seiradakis et al. 1989). Assuming a GC distance of 
8.3kpc, 0.5°corresponds to a projected distance of ~ 70pc. Even PSR J1745—2900, which overlaps the 
position of Sgr A* on this scale (~ 3"= 0.1 pc offset), is too distant for gravity tests. Figure courtesy of 
B. Klein. 


reactions should occur as they do in globular clusters. Voss & Gilfanov (2007) show that tidal 
capture of NSs by low-mass main-sequence stars is the dominant mechanism for forming binary 
NS in approximately the central l 7 of M31 and these can produce millisecond pulsars (MSPs) 
via long-term accretion. Along with similarly captured stellar-mass BHs, accreting objects in the 
central part of M31 can account for the excess number of point X-ray sources in the inner bulge of 
M31. The same processes should occur in the GC over a similar-sized region (0.1-0.2 kpc) (Muno 
et al. 2005). 

Arguably the strongest evidence for a significant NS population in the central parsec now 
comes from the recent detection of a radio loud magnetar - PSR J1745—2900, an inherently rare 
class of pulsar, just ~3"(0.1 pc) from Sgr A* (Kennea et al. 2013; Mori et al. 2013; Eatough 
et al. 2013a; Shannon & Johnston 2013). Along with the five pulsars previously discovered in 
proximity to Sgr A* (at 10 7 —15 7 - see Figure |2| (Johnston et al. 2006; Deneva et al. 2009)) these 
objects cannot be explained as part of a Galactic disk population because the number of foreground 
pulsars expected in the pulsar surveys to date is much less than one. Consequently, the question of 
the GC pulsar population has been recently revisited by a number of authors (Chennamangalam & 
Lorimer 2014; Zhang et al. 2014; Dexter & O’Leary 2014). As will be described in Section ||, it is 
pulsars in compact orbits (P or b < 1 yr - corresponding to a semi-major axis of ~ 160 au) around 
the BIT that will enable the most precise tests of GR. Recent simulations suggest that up ~ 200 
pulsars could orbit within 4000 au of the BH, the closest of which could have a semi-major axis of 
~ 120 au (Zhang et al. 2014). It is only through observations with the SKA that the full GC pulsar 
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population (whatever its size) will be detected. 

Certain dark matter models predict that dark matter particles can efficiently accumulate inside 
compact objects such as NS, eventually causing them to collapse into a BH. The GC pulsar popu¬ 
lation, as detected by the SKA, can thus also be used to place stringent constraints on dark matter 
models, with older pulsars and MSPs being the most constraining (Bramante & Linden 2014). 

3. Testing gravity with GC pulsars 

With their large moments of inertia and fast spin rates, radio pulsars are excellent clocks. The 
clocklike property of pulsars orbiting Sgr A* will provide unique probes of the spacetime, plasma 
and dark matter around the massive central BH and are a powerful probe of GR (Paczynski & 
Trimble 1979; Wex et al. 1996; Pfahl & Loeb 2004) itself. As we show below, GR and other 
perturbations of GC pulsars are large enough so that even rather “noisy” pulsars, relative to typical 
precision timing requirements, will do the job. 

Pulse times-of-arrival (TOAs) from an orbiting pulsar can be described in terms of the stan¬ 
dard Keplerian parameters and in a model-independent parameterized post-Newtonian (PPN) for¬ 
malism. Any deviations from the PPN values expected from GR would provide strong evidence 
of a post-GR theory of gravity. The most recent illustration of this method of testing theories 
of gravity is from the double pulsar PSR J0737—3039A/B, for which the measured PPN param¬ 
eters indicate agreement with GR at better than the 1% level (Kramer & Stairs 2008). How¬ 
ever, PSR J0737—3039A/B still represents a relatively weak field case. Consider the depth of 
the gravitational potential to be GM/ac 2 , for a total system mass M and semi-major axis a, with G 
and c being the usual physical constants. The current population of binary pulsars probes only to 
GM/ac 2 < 10 5 5 . By contrast, a pulsar with an orbit comparable to the star SO-2 (~ 15 yr) would 
probe to GM/ac 2 ~ 10 45 ; pulsars with shorter orbital periods could probe to even larger values 
of GM/ac 2 (Figure |]). Star-star interactions of course can cause perturbations that make GR tests 
more difficult. 

For a pulsar in a close orbit (P ol -b < 1 yr) around Sgr A* one expects to see a number of 
relativistic effects, like time dilation, Shapiro delay and relativistic peri-center precession, all of 
them with a magnitude, that as we mention above, does not even require a particularly high timing 
precision to measure them with excellent accuracy (Liu et al. 2012). The time dilation, for instance, 
is expected to be of order 10 to 20 minutes. Consequently, soon after the discovery relativistic 
effects will be well measured, give the mass of Sgr A* with high precision, and allow for the 
first gravity tests. After a few orbits, one will be able to separate the Lense-Thirring drag from 
the orbital motion and measure the spin of the BH with a fractional precision of 10 2 to 10 1 
(Liu et al. 2012). A first test is to see if the dimensionless spin-parameter does not exceed one, 
as expected for a Ken - BH in GR (cosmic censorship conjecture - see Liu et al. (2012)). The 
quadrupole moment of Sgr A* leads to a characteristic signature in the timing residuals, directly 
related to the quadrupolar structure of its potential (Figure |). By this, the effect of the quadrupole 
can uniquely be identified in the motion of the pulsar, and its magnitude can be measured with 
good to very good precision, allowing for a test of the no-hair theorem due to the unique relation 
of mass, spin and quadrupole for a Kerr solution (Wex & Kopeikin 1999; Liu et al. 2012). Figure |] 
illustrates some of the relativistic effects that could be determined from measuring pulse TOAs of 
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time (yr) 

Figure 3: Characteristic pulsar timing residuals caused by the quadrupole moment of the Sgr A* BH in 
simulated times-of-arrival (TOAs) over three orbital cycles (Liu et al. 2012). The orbital period is assumed 
to be /(,,-h = 0.3yr and the eccentricity e = 0.5. The red line is the fitted pulsar timing model. 

a pulsar orbiting Sgr A* in extremely short orbits ( P m ^ < 100 hr). For instance, measurement of 
the relativistic apsidal advance, the second order Doppler shift, and the Shapiro delay (not shown) 
yield redundant constraints on the individual masses in the system as well as provide precision tests 
of GR. The Shapiro delay is detectable for elliptical orbits even when the orbit is face on. Geodetic 
precession of the spin axis will cause secular or periodic changes in measured pulse shapes if 
there is misalignment of the spin and orbital angular momenta, and secular orbital decay from 
gravitational radiation will be measurable for the more compact orbits (F or b <10 hr). Moreover, 
as Figure [| demonstrates, only for these shortest orbital periods will the orbital decay time become 
so small that we will have little likelihood of detecting any such pulsar. Perhaps the most notable 
aspect of GC pulsars is that orbital perturbations are large enough so that even canonical (10 12 G) 
pulsars with long periods (> 0.5 s), which ordinarily are considered to be noisy clocks, are suitable 
for making timing measurements. This means that any pulsar found in the GC is useful for the 
SKA pulsar program. 

While many effects contribute to imprecision of TOAs (Cordes & Shannon 2010), the dom¬ 
inant ones for GC pulsars with canonical periods will be intrinsic spin noise (Hobbs et al. 2010; 
Shannon & Cordes 2010) and the large pulse widths that result from scattering broadening. With 
a conservative goal of about 1 s timing “precision” (in quotes because detection of gravitational 
waves with millisecond pulsars requires < 100 ns precision - Janssen et al. (2014)), many GR 
effects and orbital perturbations from stellar encounters can be measured one to two orders of mag¬ 
nitude better than the best IR interferometric measurements. The RMS spin noise for most objects 
will be below the 1 s target and the timing error from the pulsar width W is 8t ~ W/(S/N), where 
S/N is the peak signal to noise ratio of the pulse. For widths W < 1 s, it is reasonable to expect 
8t ~ 0.01 to 0.1 s, providing even more precise orbital measurements. 

At present, all these considerations rest on the assumption that the environment around Sgr 
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Figure 4: Relativistic effects that 
could be determined from measuring 
pulse TOAs from a pulsar closely or¬ 
biting Sgr A*. The bottom axis shows 
the orbital period, and the top axis 
shows the semi-major axis of the or¬ 
bit, in units of the Schwarzschild ra¬ 
dius of Sgr A*. All of the effects 
shown, except for frame dragging, 
have been detected in timing observa¬ 
tions of binary pulsars. In all panels, 
the solid line indicates a circular orbit 
(e = 0) while the dashed line indicates 
a highly eccentric orbit (e = 0.9). For 
reference, the apsidal advance rates 
of PSR J0737—3039, B1913+16, and 
Mercury are too small to be plotted 
(Cb < 20degyr *). 


A* is sufficiently clean, i.e. external perturbations on the pulsar motion are negligible. As shown 
in Figure [j] the matter content of the GC consists not only of Sgr A* , but early-type, late-type 
and evolved stars, dark matter, and magnetized plasma all of which will cause many perturbations 
that can be detected in pulsar TOAs. By measuring the Newtonian precession of elliptical orbits, 
estimation of the smoothly distributed mass, including both stars and dark matter, will be possible 
(Weinberg et al. 2005). In addition, strong constraints can be made on the density of stars from two- 
body perturbations of pulsar orbits (Weinberg et al. 2005; Merritt et al. 2010). To model potential 
external perturbations to TOAs, and see how they can be accounted for, or how much they degrade 
our capability to test GR, is subject to current research. 

Chromatic effects from intervening magnetized plasmas are discussed in the next section; these 
effects are strongly frequency dependent and can be removed (after extracting all valuable scientific 
information on the properties of the GC interstellar medium) through appropriate multi-frequency 
measurements. 

4. Chromatic radio propagation effects 

Individual pulsar TOAs will show the well known dispersion and multipath delays caused by 
the integrated electron column density (the dispersion measure, DM) and by scattering from small- 
scale (< 1 au) density irregularities, characterized by the scattering measure (SM). Polarization 
measurements of pulse profiles will show Faraday rotation from which the rotation measure (RM) 
can be determined, thus constraining the integrated magnetic field. 

Chromatic perturbations from magnetized plasma along the line of sight include the differen¬ 
tial dispersion delay across a narrow bandwidth Av (in MHz), Af(v) = 8.3 jtsDMAv/v 3 , where v 
is the observing frequency (in GHz) and DM is the dispersion measure (in pc cm -3 ), which is the 
electron column density, DM = f dln e . Such delays are removed from multi-frequency data both 
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as part of the survey analysis (typically by the use of filterbanks) and also, with greater precision, 
during timing measurements. 

Towards the GC, a more important factor is pulse broadening from scattering caused by fluc¬ 
tuations in electron density 8n e that extend down to scales of a few hundred kilometers and have a 
roughly Kolmogorov wavenumber spectrum (Armstrong et al. 1995). For the strong scattering to¬ 
wards Sgr A*, first identified by angular diameter measurements of OF1/IR masers and background 
AGNs, the pulse broadening has a time scale z s <=< f)v" 4 SM, where D is the pulsar distance and 
the scattering measure SM is the line of sight integral of C 2 , the coefficient of the wavenumber 
spectrum for 8n e (Cordes & Lazio 2003). Unlike dispersion, which is easily removed to maximize 
the S/N of a pulse, the deleterious effect of scatter broadening cannot be corrected for when many 
“scintles 1 ” are summed across the observing bandwidth, as they will be for GC pulsars. 

Magnetized plasma Faraday rotates the linear polarization by an angle v 2 RM and is rou¬ 
tinely measured for most pulsars using full Stokes parameter pulse shape measurements. In the GC 
this effect was demonstrated by measurements of the pulse polarization of the recently discovered 
magnetar, PSR J1745—2900, that revealed the second largest RM known in the Galaxy (the largest 
being that of Sgr A* itself) (Eatough et al. 2013a; Shannon & Johnston 2013). From X-ray obser¬ 
vations of the hot gas phase that Sgr A* is accreting from (and which the magnetar is likely located 
in) a simple density profile was used to place limits on the strength of the magnetic field in this 
region; a value which can affect the dynamics of the accretion process onto the supermassive BH 
(Eatough et al. 2013a). 

Once found, pulsars in the GC will therefore provide individual values of all three measures 
(DM, SM, and RM) for different lines of sight into the GC. Though not measured directly, the 
emission measure EM (the integral of n 2 e ) can be related to DM and SM so as to constrain the 
“outer scale” of the wavenumber spectrum of 8n e (Cordes & Lazio 2002). The EM, while large 
toward Sgr A* , does not cause significant free-free absorption that would prohibit detection of 
pulsars as is clear since Sgr A* is seen through optically thin plasma down to at least 1 GHz. 

5. Finding and timing pulsars in the GC with the SKA 

5.1 Past searches and updates on GC scattering 

Despite predictions that a large number of radio loud pulsars should be present in the GC (see 
Section 2.), at the time of writing just six are known within 0.5° (Figure j|). It has been understood 
that the deficit in GC pulsars could be primarily explained by extreme scattering of radio waves in 
the GC (Cordes & Lazio 1997; Lazio & Cordes 1998; Cordes & Lazio 2002, 2003). As mentioned 
in Section 4. scattering, which causes temporal broadening of pulses and a corresponding reduction 
in pulse S/N, has a strong frequency dependence (°< U 4 ). A clear mitigation strategy for dealing 
with pulse broadening due to scattering is to observe at frequencies much higher than typically used 
for pulsar searching. The significant constraint on this mitigation strategy is that typical pulsars 
have power-law spectra ( S v <=< v a ), usually with very steep spectra indices ((a) ~ —1.7) so their 

1 Scintles are the constructive interference patches seen in the frequency-time plane for nearby pulsars. The scat¬ 
tering towards Sgr A* is strong enough to prohibit measurement of individual scintles because they are narrower than 
typical frequency channel widths even at higher radio frequencies e.g. 14 GHz. 
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flux densities decline rapidly with increasing frequency (Kramer et al. 1998). Fortunately, there is 
a subset of pulsars with relatively flat spectra. For example, of 293 pulsars with derived spectral 
indices in the ATNF catalog 2 (Manchester et al. 2005), approximately 40 (or 14%) have a > — 1. 
Moreover, a small number of pulsars have been detected at frequencies as high as 90 GHz (Kramer 
et al. 1997; Morris et al. 1997). With such pulsars in mind, surveys of the GC at progressively 
higher radio frequencies have been performed (Klein et al. 2004; Johnston et al. 2006; Deneva 
et al. 2009; Deneva 2010; Macquart et al. 2010; Bates et al. 2011; Eatough et al. 2013b; Siemion 
et al. 2013) and resulted in the detection of five of the six pulsars displayed in Figure [2] (Johnston 
et al. 2006; Deneva et al. 2009); albeit at lower frequencies < 3 GHz. 

The sixth known GC pulsar, PSR J1745—2900, was first identified at X-ray wavelengths (Ken- 
nea et al. 2013; Mori et al. 2013) before detection in the radio (Eatough et al. 2013a; Shannon & 
Johnston 2013). Recent measurements of the pulse broadening in this pulsar indicate the degree of 
scattering toward the GC is less than previously predicted (Spitler et al. 2014). In addition, mea¬ 
surements with the Very Long Baseline Array (VLBA) of the angular scatter broadening show that 
it has an angular size equivalent to Sgr A* itself. This suggests that PSR J1745—2900 and Sgr A* 
are behind the same scattering region. This fact coupled with the smaller pulse broadening time 
indicate that the scattering region (if modelled as a single thin screen of material) could be distant 
from the GC (Bower et al. 2014). However, the authors note that the paucity of GC pulsars still 
indicates that regions of intense scattering might exist in the GC. These regions could be of a more 
complex structure than in the simple picture described above (Spitler et al. 2014). Therefore high 
observing frequencies (> 9 GHz) might still be necessary to detect GC pulsars. 

In the following section we use the latest information about GC scattering to re-investigate the 
prospects for GC pulsar searches with the SKA. 


5.2 Prospects for GC pulsar searches with the SKA 

Our top goal is to find and then time radio pulsars within the central parsec (30" radius), 
with the highest priority placed upon pulsars as close as possible to Sgr A* . A secondary goal is 
to conduct a more complete census of the larger region (~ 30') in order to understand better the 
overall GC NS population. While the pulse broadening time due to scattering appears to be less 
than previously predicted, the effect is still large enough to necessitate pulsar searches of the GC at 
higher frequencies. In Figure [5] we show the detection sensitivity of pulsar periodicity searches of 
the GC performed with SKA 1-MID and utilizing three of the proposed frequency bands outlined 
in (Dewdney et al. 2013). The figure shows the radio “pseudo luminosity” at 1.4 GHz (given by 
L p = .S' 14 ooD 2 . where .S'i 400 and D are the period-averaged flux density at 1.4 GHz and the distance 
respectively) versus the spin period. The dots are the actual values for known pulsars in the current 
(May 2014) ATNF catalog (Manchester et al. 2005). The curves plotted for each frequency band 
are derived from estimations of the minimum detectable flux density, .S mm , given by the radiometer 
equation as applied to pulsar observations, 


c _ a S/N m j n 7 sys I W e ff 

min “ 1 GV'AW Av T obs V flk’ 

2 http://www.atnf.csiro.au/people/pulsar/psrcat/ 


(5.1) 


9 








Galactic Centre Pulsars with the SKA 


R. R Eatough 


where /3 is a digitization factor, S/N m j n is the detection threshold, r sys is the system temperature, 
G is the telescope gain (in units of KJy '), N po i is the number of polarization channels summed, 
Av is the bandwidth, T a b s is the observation time, P is the pulsar spin period and W e ff is the ef¬ 
fective pulse width. Certain parameters have fairly standard values, namely, /3 « 1, /V po | = 2, and 
S/N min ~ 10. The effective pulse width is given by the vector sum of the intrinsic pulse width and 
any unaccounted for pulse broadening effects in the following: 

Wrff = y f (Wg 2 0 + T^ amp + T5 m + T s 2 ) (5.2) 

where W 50 is the intrinsic pulse width at half the pulse height (here assumed to be 0.05 P), T samp 
is the data sampling interval, Tdm is the intra-channel dispersion broadening (in a filterbank) and 
T s is the scatter broadening. The curves, which mark the minimum detectable pulsar spin period 
and luminosity, are plotted for pulsars at the distance of the GC and take into account the effects 
of scatter broadening, based on measurements of T s from PSR J1745—2900 (Spider et al. 2014), 
and assume a typical spectral index of a = —1.7, thereby allowing the limits to be scaled to the 
equivalent sensitivity at 1.4 GHz. Further details of the numbers used in defining these curves are 
given in the caption to Figure |5| 

If the scattering toward PSR J1745—2900 is representative of the inner GC as a whole Figure || 
shows that the common pulsars, with spin periods around 0.5 s, might be detectable at frequencies 
as low as 2.4 GHz (Band 3) with SKA1-MID. For mildly recycled pulsars, with spin periods from 
a few milliseconds to a few tens of milliseconds, 4 GHz (Band 4) observations would be needed. 
Frequencies up to at least 9 GHz (Band 5) might be required to detect the longer lived MSPs in the 
GC. As well as the possibility of detecting MSPs, GC pulsar surveys at higher frequencies offer the 
additional benefit of reduced background sky temperatures, improving sensitivity. 

For comparison with current facilities a line for a 9 GHz GC search with the Green Bank Tele¬ 
scope (GBT) has also been marked. The fact that a large fraction of the known pulsar population 
is already above this line might be of cause for concern for the prospects for finding more GC 
pulsars, however, (Chennamangalam & Lorimer 2014) note that current non-detections from past 
GC surveys cannot yet rule out a large GC pulsar population which still need deeper surveys to be 
detected. SKA1-MID and the SKA are needed to address this problem. 

Our search methodology consists of both periodicity searches, including correction for orbital 
motion, and searches for transient sources, because some pulsars are more easily detectable in 
single-pulse detection algorithms than in Fourier searches. In addition, the potential for discovering 
new classes of objects is maximized by searching for individual bursts. 

Periodicity Searches: Periodicity searches arc a standard tool for searching for pulsars and have 
been used for many years, in many directions, not limited to the GC. These algorithms 
include standard de-dispersion to remove dispersion delays, orbital demodulation, Fourier 
analysis of the resulting time series, and identification of signals above threshold, which 
requires classification as interference, instrumental contamination, or pulsars. For the GC, 
the most important aspect of algorithmic development concerns orbital demodulation, as it 
is for the shortest orbital period pulsars that the maximum GR effects can be measured. The 
simplest orbital demodulation is an acceleration search, which removes a parabolic signature 
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Figure 5: Detectability of pulsars in the GC via a periodicity search of a 6 hr integration with the GBT at 
9 GHz and SKA1-MID at three of the proposed frequency bands. Points denote the period-averaged “pseudo 
luminosity” at 1.4 GHz (L = S wo&D 2 ) versus spin period for known pulsars in the ATNF catalog (Manchester 
et al. 2005). The curves show sensitivity limits derived from the minimum detectable flux density (Equa¬ 
tion 5.1). Pulsars are detectable if they are above a given curve. We have used a GC distance of 8.3 kpc; 
a S/N = 10 detection threshold; a bandwidth Av = 0.8 GHz for the GBT, and bandwidths Av = 1.4, 2.4 
and 9.2 GHz for SKA1-MID Bands 3, 4 and 5 respectively (Dewdney et al. 2013); a system temperature, 
neglecting the effects background sky temperature, of T s ys =27 K; and a spectral index of a = —1.7, while 
recognizing that some objects have flatter and some steeper spectra. Pulse broadening, z s , values in the 
direction of Sgr A* are provided by measurements from the GG magnetar PSR J1745—2900 (Spider et al. 
2014). The curves will move lower for pulsars having spectra less steep than v _11 . The curves will also 
move lower if multiple integrations are combined (e.g., through incoherent combination of power spectra 
(Eatough et al. 2013b)). On the other hand, if for example the sensitivity of the core of SKA1-MID that 
is used for the pulsar search is reduced, by making the collecting area two thirds of the proposed size, the 
curves will move upward by an amount of log(1.5). 


as an approximation to a Keplerian orbit. This approach works for data sets that have a 
duration which is a small fraction of the orbital period (Ransom et al. 2003). 

Transient Searches: For compact orbits around Sgr A* , such as those < 100 hr shown in Figure |j. 
the large GR effects that make such pulsars extremely attractive targets will also make them 
more difficult to detect. Fast geodetic precession, for example, may mimic an orbital period 
over a 6 hr observation time and may be detectable at the output of an acceleration search. 
However, precession of the pulsar beam will make objects intermittent because the beam 
will precess out of our line of sight for a sizable fraction of the precession period. This 
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effect has been seen in several of the NS-NS binaries, including the Hulse-Taylor pulsar 
(Kramer 1998). Consequently, single epoch search observations may be inadequate to find 
such sources. Moreover, some pulsars emit “giant” pulses, with strengths that arc 1000 x 
the mean pulse intensity (Hankins & Rickett 1975; Hankins et al. 2003; Cognard et al. 1996; 
Johnston & Romani 2003). In the case of the Crab pulsar, some giant pulses outshine the 
entire Crab Nebula (Hankins & Rickett 1975), corresponding to brightness temperatures of 
order 10 31 K; (Hankins et al. 2003) observed “nano-giant” pulses having durations of only 
2 ns and implied brightness temperatures of 10 38 K, by far the most luminous emission from 
any astronomical object. 

Once GC pulsars are found with the SKA, our approach will be similar to that already used 
by e.g. Deneva et al. (2009), namely they will be timed over short and moderate time periods 
to characterize their spin properties and any orbital motion. Each will be observed at multiple 
frequencies to find the optimal band for timing precision, which depends on period, pulse width, 
flux density and spectral index, and on the exact level of scattering. Depending on the object, we 
will then embark on a long-term timing campaign for measuring GR effects and probing the GC 
environment, with a cadence determined by the orbital period around Sgr A* and on any stellar 
orbital motion. 
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